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Abstract 

A perturbation approach is used for analysis of a near-cloak in shielding a finite scatterer 
from an incident flexural wave. The effect of the boundary conditions on the interior surface 
of the cloaking layer is analysed in detail, based on the explicit analytical solutions of a wave 
propagation problem for a membrane as well as a Kirchhoff flexural plate. It is shown that 
the Dirichlet boundary condition on the interior contour of the cloak significantly reduces 
the cloaking action in the membrane case, and it also makes cloaking impossible for flexural 
waves in a Kirchhoff plate. 


1 Introduction 


The cloaking of acoustic and electromagnetic waves has been extensively developed in terms 
of theoretical design and practical implementation in papers [6}[l5||l6,19, 23][27, 28f30j. These 
publications are based on the transformation optics approach, that may employ a non-conformal 
push-out map. The theory equally applies to linear waterwave problems, as discussed in [ll . 

Elastic cloaking, both for vector problems of elasticity and for flexural waves in Kirchhoff 
plates, brings new challenges regarding the physical interpretation of equations and boundary 
conditions in the cloaking region. Indeed, this is now well understood in the context of a lack 
of invariance of the governing equations of elasticity with respect to a non-conformal cloaking 
transform. Analysis of elastic cloaking problems had been reported in 3, 22j[2lj[25] . 

In the recent papers [Ij[7j 12j 1.3. 2?)]. an advance has been made in the theoretical analysis, 
design and the physical interpretation of a cloak for flexural waves in Kirchhoff plates. Although it 
is known that membrane waves (solutions of the Helmholtz equation) drive the wave propagation 
in Kirchhoff plates (see, for example, [T 2fop()[2T ), the interface and boundary conditions provide 


coupling between solutions of the Helmholtz and modified Helmholtz equations in problems of 
scattering of flexural waves. In [4j|7| , we have analysed in detail the transformed plate equation, 
and have shown that in the cloaking region the physical interpretation requires the presence of 
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prestress and body force terms. Subject to such an allowance, the cloaking of defects in Kirchhoff 
plates becomes feasible and well understood. In the present paper, we make an emphasis on the 
important issue of the choice of boundary conditions, which are set at the interior contour of the 
cloak. 

In cloaking transformation problems, it is common that boundary conditions on the inte¬ 
rior boundary of a cloaking region are not addressed 12 J13 . This matter is rarely discussed 


and numerical simulations for cloaking are commonly presented without additional comments 
regarding these boundary conditions. The reason is that for a singular map, which “stretches” 
a hole of zero radius into a finite disk, in the unperturbed configuration there is no boundary. 
Nevertheless, in the subsequent numerical computations the singularity of the material constants 
in the cloaking region is replaced by regularised finite values. Namely, the boundary conditions 
are required for the numerical computations and usually natural boundary conditions that follow 
from the variational formulation are chosen, i.e. these are Neumann boundary conditions on the 


interior contour of the cloaking region 24 26 


One can ask a naive question related to an illusion rather than cloaking. For example, if 
one has a carrot cake, would it be possible to make it look like a fairy cake instead. In turn, 
could one make a large void in a solid look smaller? Purely naively, the latter can be addressed 
through a geometrical transformation in polar coordinates: 

r = oq + a 2 R, 9 = 0, 


where 


and 


ai = 


I?2(-Rl — O’) 

R ‘> - a 


OL 2 = 


I?2 — Rl 
Rn — a 


a < R < R 2 , R\ < r < R 2 - 


with a being a small positive number and ( R, 0) being the coordinates before the transformation 
and (r, 9) after the transformation, and R\,R 2 being positive constants, such that R\ < R 2 . 
Such a transformation, within the ring Ri < r < R 2 , would correspond to a radial non-uniform 
stretch, and assuming that outside the disk r > R 2 there is no deformation (i.e. r = R), we 
obtain a coating that would lead to an illusion regarding the size of a defect in a solid. Push-out 
transformations of such a type have been widely used, and have been applied to problems of 


cloaking, for example, in papers 14 19 


Formally, if u represents an undistorted field in the exterior of a small void, of radius a and 
we assume the series representation 


OO 

u{R,0) = Y V n (R)W n (Q), 

n =—00 


with V n and W n being basis functions, then the distortion, represented by the fields u\ and U 2 
can be described through a “shifted” series representation, as follows 

OO 

u = ui(r,9)= Y v n 

n =— 00 

and 

OO 

u = U 2 {R, 0 ) = 'Y / Vn(R)W n ( 0 ), when r > R 2 - 

n =— 00 


r — ol\ 
Ot-2 


W n {9), when Ri < r < R 2 
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Such a transformation delivers an illusion, which makes a finite circular void of radius R±, look 
like a small void of radius a. 

The above argument may appear to be simplistic, but we are going to show that it works 
exactly as described, for a class of problems governed by the Helmholtz operator and by the 
equations of vibrating Kirchhoff plates. In this cases, the basis functions V n and W n are chosen 
accordingly, and are written in the closed form in the main text of the paper. 

We consider the problem in the framework of singular perturbations and, instead of regu¬ 
larising the singular values of material parameters after the cloaking transformation, we begin 
by introducing a small hole and apply the cloaking transformation to a region containing such 
a small hole with appropriate boundary conditions already chosen. This is consistent with the 
analysis presented in [8 

In the present paper, we show that the cloaking problem is closely related to that for an 
infinite body containing a small hole of radius a where we set a boundary condition of either 
the Dirichlet (clamped boundaries) or Neumann (free-edge boundaries) type. This is done for 
both membrane waves governed by the Helmholtz equation and flexural waves that occur in 
Kirchhoff plates. Analytical solutions are presented for the cloaking problems and it is shown 
that the degree of cloaking is highly dependent on the boundary condition on the interior contour 
of the cloaking region. In a plate, in the presence of an incident plane wave along the a;-axis, 
the scattered flexural field, u s , outside the cloaking region, when a —> 0 with a clamped interior 
boundary, has the asymptotic representation at a sufficiently large distance R from the centre of 
the scatterer 
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a) 


where /3 4 = phu 2 / Dq with radian frequency w, plate flexural rigidity Do, plate thickness h and 
density p. This immediately suggests the absence of any cloaking action as the above asymptotic 
representation corresponds to a finite point force initiated by a rigid pin at the origin (see EH)- 


On the contrary, the free-edge interior boundary in the cloaking layer for a flexural plate produces 
high-quality cloaking action, and the corresponding asymptotic representation of the scattered 
field u s at sufficiently large (3R becomes 


u s ~ \J exp (i (fiR - 0) 4 (^™^) W a ) 2 = °{P a Yi as P a -> 0 and /3R -> oo. (2) 

This confirms that the regularisation algorithm that refers to a small free-edge hole produces a 
scattered field proportional to the area of this small hole, and it tends to zero as /3a —► 0. 

For the “cloaking type” map considered in this paper, the original hole does not have to be 
small, and we will refer to a “cloaking illusion”, which results in an object being mimicked by 
an obstacle of a different size. 

The structure of the paper is as follows. In Section [5] we introduce the notion of a cloaking 
transformation. Section [3] addresses the singular perturbation approach for an elastic membrane. 
In particular, in Section [3.1| we include an analytical solution, together with asymptotic estimates, 
for the model problem of scattering of membrane waves from a small circular scatterer. Section 


the full cloaked problem. Section [4] presents the analytical solution and the asymptotic analysis 
of scattering of flexural waves in a Kirchhoff plate for the biharmonic cloaking problem. The 
concluding remarks are included in Section [5j where we summarise the fundings regarding the 
cloaking action for different types of boundary conditions at the interior boundary of the cloaking 
region. 


3.2 shows the relationship between the model problem for a body with a small scatterer and 
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2 Cloaking transformation 

We aim to consider flexural waves in thin elastic plates. As shown in @ 0 . after the cloaking 
transformation such waves can be interpreted as time-harmonic flexural displacements in a pre¬ 
stressed anisotropic thin plate. Firstly, we review the cloaking transformation procedure and 
then we discuss the transformed equations and their physical interpretation. 


2.1 Push-out transformation 


A non-conformal transformation is introduced to define a cloaking coating. We use the notations 
X = ( R , 0) T and x = (r, 9) T for coordinates before the transformation and after the transfor¬ 
mation, respectively. Here ( R , 0) are polar coordinates. With reference to 14 19 23 27 28 , we 
use the radial invertible “push-out” map T : X —> x defining the new stretched coordinates 
(r,8) as follows. When |X| < R 2 , the transformation x = T(X) is given by 


r = Ri+ (fl2 p r ^ r, 9 = e, (3) 

so that i?i < r < i ?2 when 0 < R < R 2 , with Ri and R 2 being the interior and exterior radii of 
the cloaking ring, respectively. 

In the exterior of the cloaking region, when JA > R 2 , the transformation map is defined as 
the identity, so that x = X. 

In the cloaking region, the Jacobi matrix F = DT/DX has the form 


R 2 — R\ 

F = - e r 

i?2 


i?2 — R\ 


R 2 r - Ri 


eg <g> eg, 


( 4 ) 


where e r = ep>, eg = e© is the cylindrical orthonormal basis and g stands for the dyadic product. 

After the above transformation, the classical Laplace’s operator will be modified accordingly. 
For a given real a, which is identified with i?i for transformation ([3]), a new differential operator 
is defined as 


V 2 = 


1 d 
r — a dr 


(r -“% 


d 2 


(r — a) 2 d6 2 


( 5 ) 


This is consistent with notations in the paper [4]. The operator V 2 is referred to as the ‘shifted 
Laplace operator’. As in [4], we also use the terms ‘shifted Helmholtz’ and ‘shifted modified 
Helmholtz’ for the operators V 2 + (3 2 and V 2 — / 3 2 , respectively. 


2.2 Physical interpretation of transformation cloaking for a membrane 


With the reference to the work of Norris 23 , we give an outline of the transformed equation. 
Firstly, the governing equation for a time-harmonic out-of-plane displacement u of an elastic 
membrane has the form 


(V x • mV x + pw 2 ) n(X) = 0, leK 2 , (6) 

where p, p , and to represent the stiffness matrix, the mass density and the radian frequency 
respectively. 

On application of the transformation x = X(X) within the cloaking region, the transformed 
equation becomes 

(v • pC(x)V + u{x) = 0, (7) 
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where 


F = Vx®, J = det F. 


( 8 ) 



We note that equation Q, similar to ([6]), describes a vibrating membrane, but with different 
elastic stiffness, which is inhomogeneous and orthotropic, and a non-uniform distribution of mass 
across the transformed region. 

In contrast, for the model of a flexural plate, it was shown in j4j,7] that an additional pre-stress 
and body force are required to provide a full physical interpretation of the equations of motion in 
the transformed region, which again corresponds to an inhomogeneous and orthotropic material. 
The fourth order governing equation for the plate is to be discussed in section[4j Special attention 
will be given to boundary conditions on the interior of the cloaking region. 

Prior to the geometrical transformation, there was no need to prescribe any boundary con¬ 
ditions at the origin. On the other hand, after the transformation, one has an interior boundary 
of a finite size, and boundary conditions are required. In particular, this issue always occurs if 
numerical simulations are carried out for an invisibility cloak using a finite element method. On 
many occasions, such boundary conditions are chosen as natural boundary conditions and this 
means conditions of the Neumann type, which correspond to a free boundary at the interior of 
the cloaking layer. 

The question arises as to whether or not imposing a Diriclilet condition on the interior bound¬ 
ary of the cloaking layer would make a difference in a numerical simulation or an experimental 
implementation of the cloak. The answer to this question is affirmative, and it is linked to the 
analysis of a class of singularly perturbed problems introduced below. 

3 Singular perturbation problem in a membrane 

It will be shown that a regularised cloak in a membrane mimics a small circular scatterer and 
furthermore, if the radius of such a scatterer tends to zero the perturbation of the incident 
field due to an interaction with the obstacle reduces to zero as well. However, the asymptotic 
behaviour as the radius of the scatterer a —> 0, does depend on the type of boundary conditions 
on the scatterer. 

3.1 Model problem: scattering of a plane wave by a circular obstacle 
in a membrane 

Consider an unbounded homogeneous isotropic elastic solid with shear modulus p, and density p 
containing a circular void of radius a. Further, incident time-harmonic elastic out-of-plane shear 
wave of frequency oj results in a total scattered time-harmonic field of amplitude w(X) which 
obeys the Helmholtz equation outside the scatterer: 

(V 2 + A; 2 ) u(X) = 0, (9) 

where k = w-Jp / p. 

We choose polar coordinates (R, 0) such that the direction of the incident plane waves is 
along the 0 = 0 line. The total displacement field, obeying the radiation condition at infinity, 
may be written as 


u(R,Q) 


J2 \i n Jn(kR) 


i n p n H^(kR) 


e ine 


( 10 ) 
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where J n is the Bessel function of the first kind and H. n l 1 is the first Hankel function. Use has 
also been made of the multipole expansion of the plane wave 


AkR cos 0 


e".““= 53 i n Jn{kR)e ine (11) 


n=— oo 


The coefficients p n are determined by the boundary condition on the circular scatterer. They 
are given by: 

J n {ka ) 


Pn = 


Hn\ka) 

Jn(ka) 

H ( n iy {ka) 


: if u = 0 on R = i 


du 

■At -=0 on R = a 


( 12 ) 


Thus the field may be written as 


u(R,Q) = e lkRcosQ -poH^\kR) - ]T i n H^(kR)(p n e me + P - n e~ in& ) 


(13) 


n —1 


For the Dirichlet boundary condition u = 0 on R = a, this may be re-written as 


u(R,Q) = e ikRcose - -£W H W {kR) g 2^-|^tf«(fci?)cos(n©) (14) 

Hq (ka) n=1 Hn ' (ka) 

Correspondingly, for the case of the Neumann boundary condition du/dR = 0 on R = a, we 
have 


i(R,0) = e 


ikR cos © 


Jo(ka) 

H^'(ka) 


H^\kR) - g 2 i n g (fca) H^(kR)cos(nQ) (15) 


(ka) 


In the asymptotic limit, when la« 1 the coefficients near Hq X \ kR) in the above formulae 
reduce to 

m ’ 0, (16) 


Jo(ka) 


and 


H^ika) 

Jo(ka) 

H^\ka) 


2 log(fca' 


as ka 


- (ka) as ka 

4 


0, 


(17) 


respectively. 

We note that in both above cases, the coefficients near the monopole term tend to zero 
as ka —► 0. Furthermore, the coefficients near the higher-order multipole terms also vanish 
as ka —> 0. With that said, the monopole coefficient in (16), corresponding to the Dirichlet 
problem, is of order 0(\ log(fca)| -1 ), i.e. decays logarithmically slowly. Thus for a given size of 
small scatterer, the scattered field is reduced when Neumann boundary conditions are applied 
at the edge of the scatterer in comparison with the scattered field when Dirichlet boundary 
conditions are applied. 


6 














3.2 Boundary conditions and the cloaking problem in a membrane 


We note that the components of the stiffness matrix C in ([7]) after the “cloaking transformation” 
© are singular at the interior boundary of the cloaking layer. In numerical computations, a 
regularisation is always introduced at the interior boundary to “eliminate” the singularity. 

Alternatively, one can follow the approach advocated in [8] when a ” near-cloak” is defined 
through the geometrical transformation of a plane with a small hole of radius a into a plane 
containing a ring, with the unperturbed exterior radius R 2 and an interior radius R \. Outside the 
ring, the transformation is equal to the identity. Within the cloaking region, such transformation 
is defined by 

r = a 1 + CX 2 R, 0 = 0, (18) 


where 


and 


a 1 = 


R 2 ^R\ — a) 
R 2 — a 


0-2 = 


i?2 — R\ 
R‘2 ~ « 


a < R < i? 2 , Ri < r < i? 2 - 


(19) 

( 20 ) 


Outside the ring, the governing equation is ©, whereas inside the ring the equation has the 
form |7]). For simplicity, assume that the material outside the ring is isotropic and homogeneous, 
i.e. /z = const. Then, equation ([7]) can also be written as 


[ u + = 0, 


( 21 ) 


where the operator V 2 1 is defined by (|5 1. It is also noted that aq —► R\ as a 0. 

The domain with a small hole of rSius a has a singularly perturbed boundary in the limit 
when a —> 0. Simultaneously, this is also a regularisation of the push-out transformation used to 
design an invisibility cloak. 

Thus, we have the following boundary value problem 


(^“ 1 + 5) Ul M )=0, R\ < r < R 2 (22) 

(V 2 +k 2 ) u 2 (R,Q) = 0, R>R 2 (23) 

where k = lo^J p//i. The boundary conditions considered on the interior contour are the Dirichlet 
boundary condition (clamped boundary) 


ui(Ri,Q) = 0, 


or the Neumann boundary condition (free boundary) 

^(Ri,9) = 0, 


together with transmission conditions at the interface r = R 2 

Ml(i? 2 , 9 ) = U 2 (7?2 5 0 ), 


(24) 


(25) 


(26) 
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(27) 


Oi-2 ' 


du\(R 2 ,Q) du 1 (R 2 ,9) du 2 (R 2 ,9) 


dr 


OR 


OR 


and the radiation condition as R —> oo. 


The solution of the Dirichlet problem (21), (22), (p3|, (24), (26), (27) may be written as 


Ui 


M) = E 


n ——oo 


i n Jn[ -(r -a,))- 

“2 J H^(ka 


~ H n ] ( — (r-ai)) 
I V«2 ) 


And 


and 


u 2 


(R,9)= E 


l J n (kR) - i n J ^ ka) H^(kR) 
Hk (ka) 


AnO 


(28) 


(29) 


It is important to note that the outer solution u 2 (r, 9) is the field of particular interest when 
examining the quality of cloaking. It is exactly the same as the field produced by a circular 


scatterer of radius a with Dirichlet boundary condition on the scatterer (equations (10) and 
©)• Thus the degree of disturbance of the incident plane wave is totally determined by the 
effect of this single uncloaked scatterer. In turn, this is then dependent on both the size of the 
circular scatterer and the boundary conditions on it. 


The formulae (281 and (29) give the closed form analytical solution for a regularised invisibility 


cloak subject to the incident plane wave and the Dirichlet boundary condition on the interior 
contour of the cloaking layer. 

In a similar way, when the boundary condition on the interior contour is the Neumann 


boundary condition (25), the formulae (28) and (29) are replaced by 


u 1 


M) = E 


1 7 

n 


— (r - aq) - i 
“ 2 / H, 


Ah" (A-A 

(ka) V a2 / 


Jjnika) 

(i) 


And 


and 


u 2 


(J2, ©) = e 


l J n (kR)-i n J y q) HA(kR) 


Hi 1] (ka) 


AnO 


(30) 


(31) 


The exterior field outside the cloaking layer is exactly the same as for the case of a small 
circular scatterer in section |3.1| with either Dirichlet or Neumann boundary conditions. The 
displacement fields for the homogeneous membrane with a circular scatterer and the cloaking 
problem are shown in Figu re [Tj they are given for the two types of boundary conditions. As a —> 0 
and according to (16) and (17) the quality of cloaking due to the cloaking layer with the Neumann 


boundary condition is higher than that due to the cloaking layer with the Dirichlet boundary 
condition. This is because the coefficient near the leading order in the scattered component of 
the solution is 0(\ log(fca)| -1 ) for the Dirichlet case and 0((ka ) 2 ) for the Neumann case. 


4 Singular perturbation and cloaking action for the bihar- 
monic problem 

Although it appears that both Dirichlet and Neumann boundary conditions for the membrane 
can be used for the design of the invisibility cloak, we have demonstrated that the logarithmic 
asymptotics in (16) and the power asymptotics of ( |T7| suggest that the Neumann condition 
produces a more efficient cloak compared to the Dirichlet cloak as a —> 0. 

This message is further reinforced by considering flexural waves. We will now consider the 
analytical solution to a further model problem when a plane flexural wave in a Kirchhoff plate 
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Figure 1: Membrane Problem. Displacement 77[u(.R, 0)] in the direction of the incident plane-wave 
propagation as a function of the distance R from the center of the scatterer or the cloaked region. 
Curves are given for the parameters a = 0.1 cm, Ri = 1 cm, R 2 = 2 cm, corresponding to cn = 0.947 
and a 2 = 0.526 cm, and k = 3 cm -1 , (a) Dirichlet problem: displacement in a homogeneous membrane 
with circular obstacle (grey continuous line, equation {13}) and in the cloaking problem (black dashed 
line, equations j28|29 ). (b) Neumann problem: displacement in a homogeneous membrane with circular 
obstacle (grey continuous line, equation (14) ) and in the cloaking problem (black dashed line, equations 
30|31 ). The sums in the multipole expansion representations of the displacements are truncated after 


the first 10 terms. 


is incident on a small circular scatterer together with a corresponding cloaking problem in an 
analogous manner to that of sections [3. 1 1 and [372) Below, we will consider the biharmonic operator 
instead of the Helmholtz operator. 


4.1 A model problem of scattering of a flexural wave by a circular 
scatterer 

Consider an unbounded homogeneous isotropic Kirchoff plate with flexural rigidity D 0 , thickness 
h and density p containing a clamped-edge hole of radius a. Further, incident time-harmonic 
elastic flexural plane waves of frequency u> result in a total scattered time-harmonic field of 
amplitude ■u(X), at position X. As in section 3.1 cylindrical coordinates (R,Q) will be used 


such that the direction of the incident plane waves is along the 0 = 0 line. The field u(R, 0) 
obeys the following problem outside the scatterer which is regarded as clamped, 


(A 2 - (3 4 )u(R, 0) = 0, R>a ; (32) 

Ou 

U = o, — = 0, R = a, (33) 

where /3 4 = phw 2 /D 0 . 

The problem is solved by factorising the operator into Helmholtz and modified Helmholtz 
operators in the usual way. Applying the radiation condition at infinity leads to the total field 
being composed of three contributions from the original plane wave, the Helmholtz and the 
modified Helmholtz fields. The general solution is 
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u(R, 0) = 51 \i n Jn(m+PnH^(m+QnK n (m 


in& 


(34) 


n =—oo 


where K n is the modified Bessel function. Use has also been made of the plane wave expansion in 
cylindrical coordinates given in equation Application of the clamped boundary conditions 
in equation (331 leads to the following expressions for the coefficients p n and q n 


Pn 

l 

K 'n{Pa) 

- K n (/3a ) 

-i n J n {/3a) 

Qn 

~ W[HP{Pa),K n (0a)] 

l 

1 

-1 

s 

[-i n J' n {t3a) \ 


(35) 


where 


W[HW{Pa),K n (pa)\ = H^\/3a)K' n (J3a) - H™'(J3a)K n (J3a) 


For the case of free-edge boundary, the normal component of the moment and the transverse 
force are equal to zero on the boundary of the small circular scatterer, which results in the 
following conditions to be satisfied 


d 2 u 

dR 2 


+ v 


1 du 1 d 2 u \ 
RdR + R 2 902/ 


= 0, at R = a, 


and 


where 


with 


Tr(u) = ^ 


T r (u) =0, at R = 
d f d 2 u 1 du 1 d 2 u \ 


OR \dR 2 ' RdR + R 2 d@ 2 J 

1 d 2 u 


1 8Mrt{u) 

R 90 


Mrt(u) = (1 — v) 


RdRdQ 


J_du 

iRdQ 


(36) 

(37) 

(38) 

(39) 


Application of the free boundary conditions in equation (33) instead of the clamped boundary 


conditions leads to the following expressions for the coefficients p n and q n 


Pn 

l 

- V~(K n (j3a)) -X+(K n (j3a)y 

'-i n X~{J n {pa))' 

Qn 

“ & 

X~(Hn\/3a))_ 

—i n V + (J n (/3a))_ 


where 

X ± { < g n {fia)) = —/3a(l - u)^(da) + [n 2 ( 1 - u) ± (/3a) 2 ]^ n (/3a), 
V r± C^ n (/3a)) = [±(/?a) 3 + (1 - v)n 2 /3aK(/3a) - (1 - v)n 2c tf n (l3a), 
with (pa) denoting Hn(fia) or K n (/3a), and 

^ = X-{H£Xda))V-{K n {Pa)) - V+(H^(J3a))X+(K n (fia)). 


(40) 

(41) 

(42) 

(43) 


The singular perturbation approach advocated above and the regularised push-out trans¬ 
formation work equally well for the case of a cloak re-routing flexural waves. Compared to 
the membrane waves, additional action is produced by evanescent waves; in particular, Green’s 
function for flexural waves remains bounded whereas for the membrane problem it is singular. 

We refer to the paper 18 that has addressed the classical solution of the scattering of a 
flexural plane wave. Here we summarise the results, with the emphasis on the asymptotics of 
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outgoing waves when (3a becomes small, 0 < (3a <C 1. As in 18 , the coefficient po in the 

L 

(44) 


expansion (341 of the flexural displacement has the following limit representation as /3a <C 1 

— 1, as (3a 


Po ~ —i, as pa —► 0, 
for the case of the clamped boundary (Dirichlet problem), and 


mv , 

Po ~i(TVo (,fa) 


as (3a —> 0, 


(45) 


for the case of the free-edge boundary (Neumann problem). Furthermore, the coefficients 
p n , \n\ > 1, vanish as (3a —> 0, for both types of boundary conditions considered above. We 
also note that the coefficients q n near K n terms are not significant, as these terms are exponen¬ 
tially small away from the scatterer, as (3R 1. 


For the Dirichlet problem, when the contour of the circular inclusion is clamped as in (33), 
the scattered field u K is 


«s~- 




oo. 


(46) 


For the case of a free-edge boundary condition, where both the transverse force and the moment 
vanish at R = a, the scattered field u s becomes 


7 t(3R 


exp 


(i {PR - v 


mv 


4/ / 4(1 - v) 


(/3a) 2 as j3R —► oo. 


(47) 


0, the scattered field produced 
corresponds to a finite force. In other words, for the Dirichlet boundary 


Equation (46) clearly illustrates that, even in the limit of (3a 
by the “rigid pin' 


condition, the singular perturbation of the boundary produces a finite scattered field and this 
remains finite in the neighbourhood of the scatterer even when the diameter of this scatterer 
tends to zero. On the contrary, for the case of the free-edge boundary, as shown by (47), the 
coefficient in the leading order term is 0(((3a) 2 ), and it vanishes as (3a —> 0. 

Comparison with © and (JTtJ) shows that, while for the Helmholtz operator, the monopole 
term in the scattered field always vanishes in the limit of pa —► 0 (both for Dirichlet and 
Neumann boundary conditions); in the biharmonic case, which corresponds to a Kirchhoff plate, 
the scatterer with the clamped boundary (Dirichlet boundary condition) delivers a finite non¬ 
zero point force in the limit of (3a —> 0. On the other hand, the case of a free-edge boundary 
(Neumann boundary condition) gives no scattering in the limit of (3a —>■ 0, as seen from (47). 


4.2 Boundary conditions and the cloaking problem in a Kirchhoff plate 


The problem to be solved in this sub-section is analogous to that in section 3.2 except that we 
will consider cloaking of flexural rather than out-of-plane shear waves. A “near-cloak” is again 
defined through the invertible geometrical transformation of a plate with a small hole of radius 
a, into a plate containing a ring, with unperturbed exterior radius i ?2 and an expanded interior 
radius R\ as in equations (18), (19) and 


Again, incident time-harmonic flexural waves of frequency w result in a total scattered time- 
harmonic fields of amplitude u±(r, 9) and U 2 (R, 0) inside and outside the cloak respectively. The 
material outside the ring is isotropic and homogeneous and the governing equation is as in Eq. 


(32). Inside the cloak, the governing equation is given by the transformed biharmonic equation 
and the classical Laplace’s operator is modified accordingly as in equation ([5]). The problem to 
be solved is defined as 


11 


















V 4 , - ^4 ) u i(r,9) = 0, Ri<r<R 2 


a 2 
74 o4 


(V 4 - p 4 ) u 2 (R, 0) = 0, 

For the case of the clamped interior boundary, the boundary conditions are 

Ui(Ri,d) = 0, 
dui(Ri,6) 


dr 


= 0. 


(48) 

(49) 

(50) 

(51) 


For the case of free-edge boundary, zero transverse force and normal component of the moment 
on the interior contour of the cloaking region correspond to transformed boundary operators, 
obtained from the variational formulation for the shifted biharmonic operator in (48) used within 
the cloak. They result in the following conditions to be satisfied at r = R\ 


T r = —c 


d fd 2 ui 1 du\ 1 <9 2 mi\ 

dr \ dr 2 + [r — ot\) dr + (r — au) 2 dd 2 J 


ot 2 


dWl r 


with 


and 


m rt =ai(l- v) 
'd 2 m 


1 


d 2 u\ 


1 


(r —ai) dd 
du\ 


= 0, 


(r — ai) drdO (r — ai) 2 d9 


Qio 


dr 2 


1 dui 1 9 2 wi\ 

(i— ai) dr (r — ai) 2 dd 2 ) 


= 0. 


In addition, on the interface between the cloak and the ambient material we have 

ui(R 2 ,d) = u 2 {R 2 ,Q), 
dui(R 2 ,d) dui(R 2l d) du 2 (R 2 ,<d ) 


«2' 


dr 


d 2 u 2 

dR 2 


1 du 2 1 d 2 u 2 \ 
R~dR + R 2 ae 2 ") 


dR 

d 2 u\ 

-1 

dr 2 


dR 

1 du\ 1 d 2 u \\ 

(r — a±) dr (r — au) 2 dd 2 J 


(52) 

(53) 

(54) 

(55) 

(56) 

(57) 


for the normalised normal component of the moment on the outer boundary of the cloak, and 

T r (u 2 ) = £ r (ui) (58) 

for the normalised transverse force on the outer boundary of the cloak. The transverse force 
T r (u) is defined in (38) and (39). 

We also supply the radiation condition as R —► oo, i.e. the scattered held is represented as 
an outgoing wave. 

For the homogeneous material, when a 2 = 1 and a± = 0, the above conditions imply the 
continuity of the flexural displacement and rotation, in addition to the normal moments and 
transverse forces. 

The series representation of the flexural displacement held inside the cloak and outside the 
cloaking region has the form 


Ml 


(r,0)= ^ a n J n (-^(r - a!)\ +b n H^ l 1) f-^-(r - ai)\ 

r> ——on - ' ^2, J \ ^2 J 


In 


— (r - ai) ) + d„K„ \ —{r - a{) 

Ot 2 / \Ot 2 


e me , R 1 < r < R 2 , 


(59) 
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and 


U2 (i?,0)= ]T i n J n (0R)+p n H^(0R)+q n K n (0R) e in0 , R> R 2 , (60) 


n——oo 


respectively. 

For each order n , there are six unknown coefficients. For the case of the clamped interior 
boundary r = R\ , a solution satisfying the six boundary and interface conditions (501 - (58) is 
given by a n = i n , b n = p n , c n = 0, d n = q n together with p n and q n satisfying 


Pn 

l 

K'nifia ) - K n {0a) 

-i n J n {l3a) 

_Qn_ 

W[H£\0a),K n (Pa)] 

- Hk lY (l3a ) Hn^ {0a) 

r inJ n(H 


(61) 


where 


W[HW(pa),K n (J3a)] = H^(pa)K' n ((3a) - H^\pa)K n (0a). 


Note that these equations for p n and q n are the same as in (|35j) for the case of the clamped 
boundary R = a, and the corresponding coefficients for the case of the free-edge condition are 
the same as p n and q n in equation p0| ) , respectively. 

This implies that the outer field u 2 (R, 0) for the cloaking problem is exactly the same as the 
field produced by a circular scatterer subjected to the same boundary conditions (either clamped 
or free-edge boundary) as those on the interior boundary of the cloak in the cloaking problem. 
The required solutions are given by equations (341 (43). The displacement fields are shown in 
Figure [5] 


Dirichlet 



Neumann 



(b) 


Figure 2: Plate Problem. Displacement TZ[u(R, 0)] in the direction of the incident plane-wave propaga¬ 
tion as a function of the distance R from the center of the scatterer or the cloaked region. Curves are 
given for the parameters a = 0.1 cm, Ri = 1 cm, R 2 = 2 cm, corresponding to ai = 0.947 and a 2 = 0.526 
cm, 0 = 3 cm^ 1 and v = 0.3. (a) Dirichlet problem: displacement in a homogeneous plate with circular 
obstacle (grey continuous line, equations (|34|35}) and in the cloaking problem (black dashed line, equa¬ 


tions j59|61} ). (b) Neumann problem: displacement in a homogeneous plate with circular obstacle (grey 
continuous line, equations |34|40| ) and in the cloaking problem (black dashed line). The sums in the 
multipole expansion representations of the displacements are truncated after the first 10 terms. 
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The asymptotic representation of the displacement field it 2 , for R » 1, outside the regularised 
cloak, considered here, is also the same as the one for the small inclusion with either Dirichlet or 


Neumann boundary conditions, as in section 4.1 In particular, according to the formulae (46) 


and (47), the cloaking action is non-existent if the interior contour of the cloaking layer in the 
Kirchhoff plate is clamped (Dirichlet boundary condition set on r = f?i). On the other hand, the 
cloaking appears to be efficient for the free-edge boundary, with the leading term of the scattered 
field to be of order (/3a) 2 as /3a —► 0. 


5 Concluding remarks. Cloaking illusion for the mem¬ 
brane and the Kirchhoff plate problems 

This paper dispels a common perception of an “invisibility” based on a singular geometric trans¬ 
formation used in the design of invisibility cloaks. It has been shown that the boundary conditions 
on the interior contour of the “cloaking” may influence significantly the scattered held. Singular 
perturbation analysis is the tool which enables one to see it. 

Based on the model of a near-cloak, for flexural waves in a Kirchhoff plate we have shown 
that the object surrounded by the cloak, appears as an infinitesimally small scatterer, with 
appropriate boundary conditions on its contour. The intuitive perception of a scattered held 
that vanishes when the diameter of such a scatterer tends to zero proves to be wrong for the case 
of flexural waves scattered by a clamped small inclusion. 

Even when a singular transformation is used to design an “exact” cloak, in the numerical 
implementation, or in a physical experiment, it is always regularised. As noted earlier, such 
regularisation is linked to setting boundary conditions on the interior contour of the cloaking 
region. 

For the near-cloak, the singular perturbation model shows that even for the case of an elastic 
membrane, governed by the Helmholtz equation, setting the homogenous Dirichlet boundary 
conditions would lead to a higher scattering compared to the case of homogeneous Neumann 
boundary conditions on the interior contour of the cloaking layer. This observation is fully 
supported by the asymptotic formulae for the coefficients of the leading term of the 

scattered field for the near-cloak. The coefficient p 0 near ( kr ) in the representation for the 

scattered field in the membrane, as ka —> 0, is 


P 0 

Po ' 


7TZ 


2 log (ka) ’ 



(Dirichlet problem), 
(Neumann problem). 


(62) 


Furthermore, the message regarding the choice of the Neumann boundary conditions is rein¬ 
forced by the analysis of the flexural waves re-routed by a cloak in a Kirchhoff plate. For the 
biharmonic problem describing such flexural waves, we have shown that setting the full clamping 
(homogeneous Dirichlet boundary conditions) of the interior contour would remove the cloaking 
action, even though the physical parameters of the cloak are chosen according to the push-out 
transformation ©>■ On the other hand, the homogenous Neumann boundary conditions, cor¬ 
responding to a free-edge interior contour of the cloak, would be appropriate to the cloaking 
action and hence routing membrane and flexural waves around a finite scatterer. For the Kirch¬ 
hoff plate waves, the coefficient po near Hq (fir) in the representation for the scattered field is 
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Membrane 


Plate 




Figure 3: Monopole coefficient: modulus of |p 0 |. (a) Membrane problem: |p 0 | is given for the Dirichlet 
and Neumann problems, equations {62}. (b) Plate problem: |p 0 | is given for the Dirichlet and Neumann 
problems, equations (63|. 


characterised by the following asymptotics as /3a —> 0 : 

Po ~ —1, (Dirichlet problem), 

Po ~ -r(/3a) 2 , (Neumann problem). (63) 

4(1 - v) 

The dependance of the monopole coefficient po as a function of the small parameters ka and /3a 
for the membrane and plate problems, respectively, is also shown in Figure im¬ 
practical applications of filters, polarisers and cloaks for elastic flexural waves are in the 
efficient design of the earthquake protection systems. The present paper shows that a choice of 
the type of boundary conditions on the interior contour of the cloak is essential. The results of 
the paper also have a wide range of applications concerning practical designs of cloaking shields 
for electro-magnetic and acoustic waves. 
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